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A model for ID traffic flow is developed, which is discrete in space and time. Like the cellular 
automaton model by Nagel and Schreckenberg [J. Phys. 1 France 2, 2221 (1992)], it is simple, fast, 
and can describe stop-and-go traffic. Due to its relation to the optimal velocity model by Bando 
et al. [Phys. Rev. E 51, 1035 (1995)], its instability mechanism is of deterministic nature. The 
model can be easily calibrated to empirical data and displays the experimental features of traffic 
data recently reported by Kerner and Rehborn [Phys. Rev. E 53, R1297 (1996)]. 
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Cellular automata (CA) are interesting for their speed 
and their complex dynamical behavior [ffl] , including such 
fascinating phenomena as self-organized criticality [§-|| , 
formation of spiral patterns [pi , or oscillatory and chaotic 
sequences of states [0,|§. Their enormous computation 
speed and efficiency is a consequence of the following 
properties, which are ideal preconditions for parallel com- 
puting: (i) discretization of space into identical sites x, 
(ii) a finite number of possible states f{x) (hi) the (par- 
allel) update at times T = t AT with an elementary time 
step AT, (iv) globally applied update rules, based on (v) 
short-range interactions with a finite (small) number of 
neighbouring sites. Despite these simplifications, cellular 
automata have a broad range of applications, reaching 
from realistic simulations of granular media or fluids 
(including interfacial phenomena and magnetohydro- 
dynamics), over the computation of chemical reactions 
HU, up to the modeling of avalanches j|. Their ap- 
plication to traffic dynamics jiolll has stimulated an 
enormous research activity |12-|14|], aiming at an under- 
standing and control of traffic instabilities, which are re- 
sponsible for stop-and-go traffic and congestion, both on 
'freeways' and in cities. 

Recently, Kerner and Rehborn |l5| ] have reported some 
characteristic properties of empirical highway traffic flow, 
which a realistic traffic model should display: (i) At small 
densities, traffic flow is stable, i.e., arbitrarily large dis- 
turbances of homogeneous traffic will disappear in the 
course of time, (ii) Above a certain critical density, any 
small perturbation will give rise to the formation of a 
traffic jam. (iii) Between the stable and the unstable re- 
gions, there exists a density interval beginning at about 
20 vehicles per kilometer, where traffic flow is metastable. 
That is, sufficiently small disturbances (so-called 'sub- 
critical perturbations') will fade away, whereas 'super- 
critical' perturbations exceeding a certain minimal am- 
plitude will cause a traffic jam. (iv) The outflow from 
traffic jams has a typical value which is independent of 
the initial conditions and, to a large extent, independent 
of the average surrounding traffic density. It varies only 
with road or weather conditions, and the average vehi- 
cle characteristics (regarding their lengths and acceler- 



ation capabilities), (v) The typical outflow is consider- 
ably smaller than the maximum flow and lies at about 
1600 vehicles/km for slow lanes, and on fast lanes be- 
tween 1800 vehicles/km (measured in Germany |l5)) and 
2100 vehicles/km (on Dutch motorways Jl^]). It is as- 
sociated with a density of 20 ± 5 vehicles per kilome- 
ter, (vi) Downstream jam fronts move with a velocity of 
— 15 ± 5km/h. Notice that properties (iv) and (vi) orig- 
inate from the uniform acceleration behavior of queued 
vehicles, resulting from the similar distances and veloci- 
ties that they share inside traffic jams. While the prop- 
agation velocity C of traffic jams is given by the dissolu- 
tion speed of a jam front, the outflow Q ut is related to 
the time gap T between successive departures from the 
traffic jam |l6)| . 

The cellular automaton proposed by Nagel and 
Schreckenberg |l(| meets the properties (i), (ii), and (vi), 
and most of the other properties can be reproduced by 
separate variants of it jl7] |i~9|. In particular, the continu- 
ous version by Kraufi et al. ]18| and the CA by Barlovich 
et al. |l^] seem to display metastable states. Moreover, 
the continuous version is in good agreement with em- 
pirical traffic data []l8|] . Here, we will present a "unify- 
ing" cellular automaton which, in a certain parameter 
range, reproduces all of the above properties. Remark- 
ably, the characteristic quantities can be calculated an- 
alytically. Moreover, consistent with macroscopic traffic 
models, the mechanism of traffic jam formation is deter- 
ministic Jj0||l5| ] rather than based on internal fluctuations 
( "randomization" ) |l^] . Being related to the optimal ve- 
locity model pi] , it originates from the delayed adapta- 
tion to an equilibrium velocity, which, in the instability 
region, rapidly decays with growing density [ p^pCj ] . 

We propose the following simple and fast discrete 
model that can be well calibrated to macroscopic traf- 
fic data: To maximize simulation speed, we first choose 
the time step AT of the temporal update as large as pos- 
sible. It is limited to AT w 1 s, since this corresponds 
to the safe time headway required for avoiding accidents. 
The spatial discretization AX of the road should not be 
larger than the minimal vehicle distance £. A fine veloc- 
ity discretization AV = AX/ AT is reached by taking a 



1 



D. Helbing/M. Schreckenberg: Cellular Automata Simulating Experimental Properties of Traffic Flow 



2 



fraction AX = £/n of £ (n e {1, 2, 3 . . .}). Velocity steps 
of the order 5 — lOkm/h require AX < £/2. In the fol- 
lowing, the spatial coordinate X = x AX and distances 
D = dAX are measured in units of AX, time T = t AT 
in units of AT, and any velocity in units of AV. 

Assume we want to distinguish A different vehicle 
types a € {1,2,... ,A}, e.g. cars and trucks [Q. Then, 
at each time t e {0, 1,2,.. .}, every site x G {1,2, ... ,1} 
can be in one of the states ft{x) = (a, v), where ft{x) — 
(0,0) corresponds to an empty site, and ft(x) — (a,v) 
with a > represents a vehicle of type a with velocity 
v £ {0,1,2,..., «max}- For safety reasons, the velocity v t 
must be smaller than the distance d t to the vehicle ahead. 
The states of the cellular automaton are updated in par- 
allel according to the following successive steps: First, 
each vehicle is moved by its actual velocity Vt to posi- 
tion x + Vt, which means ft+i(x + Vt) = ft(x), if a > 
and vt > 0. Then, the states of the previously occupied 
positions are reset to f t +i{x) — (0,0). Any other site 
keeps its previous state. Finally, all vehicle velocities are 
modified along with the proposed acceleration rule 



"t+i 
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J t+i 



X a [v a (d t+1 ) - w t ]j , (1) 
1 with probability p, if v' t+1 > 



otherwise, 
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where [j/J is defined by the largest integer i < y. There- 
fore, the above equation implies v t +i < X a v a(dt+i) + 
(1 — X a )v t , meaning that the new velocity is a weighted 
average of the previous velocity v t and the optimal ve- 
locity v a of vehicle type a, or somewhat less. A small 
value of the model parameter A a > relates to a great 
inertia of vehicle motion, whereas a large value A a < 1 
implies a fast adaptation to the distance-dependent op- 
timal velocity v a (d). The corresponding adaptation time 
is r a — AT/X a . If the (back-bumper-to-back-bumper) 
distance d t +\ to the next vehicle exceeds a certain finite 
value dfi n , the vehicles do not interact, and v a is given 
by the maximum velocity v^ax °f vehicle type a. For 
small distances, v a should be determined by the velocity- 
dependent safe distance d{v a ) ~ ij AX + v a , required to 
avoid accidents. The model parameter ('slowdown prob- 
ability') p describes individual velocity fluctuations due 
to delayed acceleration (imperfect driving). Here, we are 
interested in the limit p — ► 0. 

In order to compare this discrete model with the ob- 
served properties of traffic flows, it is necessary to inves- 
tigate aggregate quantities (h) . These are defined by 
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The vehicle density is given by p(x AX, t AT) = (6(a)) / 
AA, the traffic flow by Q(xAX,t AT) = (w6(a))/AT, 
and the average velocity by V(x AX,t AT) — 



Q(x AX, t AT)/p(x AX, t AT), where 6(a) = 1 for a > 
0, otherwise = 0. During the simulation runs, the den- 
sity minima and maxima, their propagation speed, and 
the associated traffic flows were evaluated automatically, 
as well as the spatially averaged vehicle velocity and traf- 
fic flow on the circular road of length L = 20 km. The 
corresponding values were averaged over several hours af- 
ter a sufficiently long transient period. Averages over the 
whole street are indicated by an overbar. 

First, let us discuss the case of one type a = 1 of vehic- 
les. Our simulation results can be summarized as follows: 
At small average densities p, homogeneous traffic flow is 
stable, and the spatially averaged velocity V is given by 
the density-independent value (w? lax — [1/Ai — 1] — p) AV. 
It is zero for high densities [Fig. 0(a)], since Ai^i (d) > 1 is 
necessary for the acceleration of a vehicle from standstill. 

At medium densities, the resulting velocity-density re- 
lation is largely dependent on Ai: Regime I: For 1 > Ai > 
A s tab(i'max>-P) (corresponding to a quasi-instantaneous 
adaptation to the optimal velocity), V{jj) is close to the 
piecewise constant relation vi(l/~p). Accordingly, the as- 
sociated average traffic flow Q is a piecewise linear func- 
tion of p. Regime II: In a certain interval A m i n (ti^ , p) < 
Ai < A m ax(Wmax;-P); traffic flow is unstable for a certain 
range p ut < 7> < Pmax of medium densities, and Q(p) 
becomes the self-organized linear relation 



T \ p jan 



(4) 



demanded by Kerner [|l6| [Figs. 0(a) and||(b)]. T denotes 
the average time gap between the acceleration of succes- 
sive vehicles. The linear relation (|]) reflects a mixture of 
free and jammed traffic with characteristic densities p ou t 
and Pjam, respectively, where the jammed regions grow 
with increasing density. Qout = Q(/Oout) is the typical 
outflow from traffic jams pq]. The slope 
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dQ 
dp 



(5) 



corresponds to their dissolution velocity |Tq| . The de- 
pendence of the spatially averaged velocity V on p is 
given by V(p) = Q(~p)/p. Regime III: In a parameter 
range A max (^ lax ,p) < A x < A sta b(t4 ax ,p), tnere is stm 
an unstable range of traffic, but Q(fi) is piecewise lin- 
ear with different slopes C and different values of pj am 
or p ou t [like in Fig. 0(b)], where the above relations are 
separately fulfilled for each linear piece. This may be un- 
derstood as crossover behavior between the cases I and 
II. Notice that the discretization of vehicle dynamics im- 
plies 1/pmax = k\ AX and T = AT/i;3 with small 
integers hi (see below). Thus, C is restricted to a few 
discrete values ^^-AV. Regime IV: For Aiw 1 ax < 1, Ai 
is so small that v\ +1 < v t which, for p ^ 0, implies that 
traffic eventually comes to rest. 
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The characteristic quantities can explicitly be calcu- 
lated. Let us show this for the simple but nontrivial 
case p —>■ and Vi(d) — min(c? — 1,3) modeling city 
traffic. With AX = 6.25 m we have Pj am = 1/AX — 
160 vehicles/km. Now, let us assume that a disturbance 
has produced a queue of vehicles with velocities v = 
and distances d = 1 to the respective vehicles ahead, and 
a free road in front of the first vehicle. We can character- 
ize the acceleration behavior of a vehicle by the sequence 

(vt',dt') (vf+udt'+i) '-^ (v t '+2,dt'+2) *-* 2 ---, (6) 

where, for A = 0.77, v t +\ — max(d t — 2,0) if (v t + 1) < 
d t < 4, and d t+1 = (d t + v* +1 - v t+1 ) being the 

velocity of the vehicle ahead) . Denoting with t' the time 
when the state (0,1) of the respective vehicle ahead has 
changed to another state, we find two alternating se- 
quences: (0,1) -> (0,1) -> (0,3) -> (1,4) -> (2,4) -> 
(2,4) -» ... and (0,1) - (0,2) -> (0,4) -» (2,4) -» 
(2, 4) — ► . . . That is, cars start to accelerate alternatingly 
every one or two time steps AT. With AT = 1 s, we find 
an average of T = 1.5 s, which implies the dissolution ve- 
locity C = — 15km/h. Moreover, the resulting density in 
front of jams is given by the evolving maximal distance 
4AX, which gives p ou t = 40 vehicles/km and, according 
to Eq. (P, an outflow of Q out = 1800 vehicles/h. All this 
is in total agreement with simulation results [Fig. 0(a)]. 

Now, let us investigate the behavior for fixed Ai = 0.77 
and density p = 80 vehicles/km, but various p [Fig. 0(a)]. 
If p is continuously reduced, the backward dissolution 
velocity C is decreasing, whereas the outflow Q ou t from 
traffic jams, and the difference between the jam density 
Pjam and the self-organized density p out downstream of 
traffic jams are rapidly growing towards an almost con- 
stant value. In particular, the amplitude [pj&m — Pout) of 
traffic jams is constant over more than five decades. 

In the deterministic case p = 0, the formation of traf- 
fic jams requires some initial inhomogeneity. Let us in- 
vestigate the response to localized perturbations of the 
form p(x, 0) = ~p + Ap{cosh~ 2 [(x — L/2)/w+] — (w + /w_) 
cosh~ 2 [(a; — L/2 — w + — W-)/w-]}, as suggested in Ref. 
p3j. Typically, one observes a piecewise linear response 
like in regime III [Fig. 0(b)]. However, its dependence on 
the perturbation amplitude Ap indicates multistability, 
i.e. the coexistence of a variety of solutions. These corre- 
spond to periodic patterns of vehicle updates (which nat- 
urally result in deterministic systems with a finite num- 
ber of states). In the presence of noise (p > 0), only one 
of them survives [Fig. |^(a)], i.e. most of them are unsta- 
ble with respect to fluctuations. This explains the role of 
randomization for the behavior resulting in regime II. 

Finally, let us focus on the density-dependent behav- 
ior for fixed Ai = 0.77 and p > 0. We find that initial 
localized perturbations of the above form, regardless of 
their amplitude Ap, are damped out for average densities 
p below some value Pci( - y ] ^ lax ,p) and above some value 



Pc4« ax ,p)- In a certain density range p C 2« ax ,.p) < 
P < Pc3(w niax ,p), the perturbation grows for any finite 
amplitude. In the density regimes Pci^maxiP) — P — 
Pc2« ax ,P) and Pc3« ax >P) < P < Pc4« ax ,p), we ob- 
serve metastability [Fig. [|, i.e. perturbations with an 
amplitude Ap > Ap cr (jj) will grow, otherwise they will 
fade away in the course of time ('local cluster effect' 
Q|3|). Notice that p c2 ~ p ou t and p c3 w p ma x. 

These findings can be understood by analogy with the 
continuous optimal velocity model by Bando, Sugiyama 
et ai. j2l| that results in the limit AT -> and AV -> 0. 
It displays stable traffic at low and high vehicle densities, 
unstable traffic on the condition dvi(d)/dd > Ai/2, and 
metastable regimes between the unstable and stable ones 
|pl| . However, in the continuous optimal velocity model, 
the jam density pj am is not independent on how a traffic 
jam is formed, since fast cars are more crowded than slow 
ones, after they had to stop. This undesired property can 
be avoided by a refined model 0] or just by a suitable 
discretization, like in the proposed cellular automaton, 
when operated in regime II. 

In summary, we have developed a cellular automaton 
for one-lane traffic which reproduces many of the em- 
piricially observed features of traffic flow in a "unified" 
way. In particular, the model showed the characteristic 
quantities of traffic flow, which we managed to calculate 
analytically. By appropriate specification of the tabu- 
lar functions v a (d) and the parameters A a , p, AT, and 
Ay, the model can be calibrated to empirical data. The 
most interesting case is to operate the model in regime II, 
since this guarantees the desired properties (iv) and (vi). 
The characteristic quantities like C and Qout are deter- 
mined by Pout, Pjam = ki AT, and T = k 2 AT/k$. The 
latter can be enforced by a suitable choice of AX and 
AT. v a (d) and A a determine the maximum velocity, the 
approximate velocity-density relation, the instability re- 
gion, and the amplitude (pj am — p ou t) of traffic jams. The 
metastable regimes and the difference between the max- 
imal possible traffic flow and the self-organized outflow 
Q ou t from traffic jams grow with increasing t)^ m . Finally, 
p allows to influence the characteristic 'wave length' be- 
tween successive traffic jams [Fig. ||(a)]. Suitable choices 
are AT G [Is, 1.3s], A a w 0.77 and p < 0.01. The opti- 
mal velocity functions v a (d) were chosen proportional to 
relations that were determined from traffic data of the 
Dutch freeway A9. The results are in good agreement 
with macroscopic traffic data [Fig. ||(b)-(d)]. 
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FIG. 2. (a) Characteristic parameters of traffic flow as 
a function of the slowdown probability p for the model 
parameters specified in Fig. hi (b) Comparison of the 
model corresponding to A„ = 0.77, p = 0.001, AT = Is, 
AX = 2.5 m, and the optimal flow relation specified by '+' 
with 1-minute-averages of single-vehicle data from the left 
lane of an undisturbed cross section of the Dutch highway 
A9 (■). Boxes illustrate the simulated average flow result- 
ing in the limit of long times, (c) Simulation of 'stop-and-go 
waves' at an average density of 28 vehicles/km for a mixture 
of 90% cars and 10% trucks (the optimal velocities of which 
are only 70% of the cars). The occuring minimal and maxi- 
mal velocities, the minimal densities, and the largely varying 
time intervals of successive breakdowns of velocity are in good 
agreement with the Dutch freeway data (1-minute-averages) 
displayed in (d). The mixture of vehicle types also explains 
the observed fluctuations in the density and average velocity 
of vehicles. 
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FIG. 1. Simulation results for vi(d) — min(d — 1,3), 
Ai = 0.77, AT = Is, AX = 6.25m, and (a) p = 0.001, 
(b) p = 0, w+ = W— = 200 m, and Ap — p/2. Illustration 

(a) shows the 'optimal flow' 7wi(l/p) (+) and the resulting 
average flow Q{p) (□) as a function of the average density ~p, 

(b) the densities inside ( — ) and in front of ( — ) traffic jams. 
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FIG. 3. (a) Critical amplitudes of localized perturbations 
with w+ — 200 m and W- = 800 m as a function of average 
density for the model specified in Fig. ^|(b). Whereas larger 
perturbations cause the formation of traffic jams, smaller ones 
will fade away (p c i = 21 vehicles/km, p C 2 = 23 vehicles/km, 
Pc3 = 150 vehicles/km, p C 4 > 164 vehicles/km), (b) The av- 
erage wave length of emerging stop-and-go waves diverges at 
p C 2 (checked for large system sizes). This is why they are not 
triggered by fluctuations below this density. 



